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ABSTRACT 

Herschel-HWl observations of high-i lines (up to y u = 10) of 12 CO, 13 CO and C 18 are presented toward three deeply embedded low-mass proto- 
stars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. 
The spectrally-resolved HIFI data are complemented by ground-based observations of lower-/ CO and isotopologue lines. The 12 CO 10-9 profiles 
are dominated by broad (FWHM 25-30 km s~') emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 
100-200 K. Several CO and 13 CO line profiles also reveal a medium-broad component (FWHM 5-10 km s~'), seen prominently in H2O lines. 
Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The 
narrow C ls O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation 
temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars. 

Key words. Astrochemistry - Stars: formation - ISM: individual objects: NGC 1333 - ISM: jets and outflows - ISM: molecules 



1. Introduction 

The earliest protostellar phase just after cloud collapse - the so- 
called Class phase - is best studied at mid-infrared and longer 
wavelengths dAndre et al.l2000h . To understand the physical and 
chemical evolution of low-mass protostars, in particular the rel- 
ative importance of radiative heating and shocks in their energy 
budget, observations are required that can separate these com- 
ponents. The advent of the Heterodyne Instrument for the Far- 
Infrared (HIFI) on Herschel opens up the possibility to obtain 
spectrally resolved data from higher- frequency lines that are sen- 
sitive to gas temperatures up to several hundred Kelvin. 

Because of its high abundance and strong lines, CO is the 
primary molecule to probe the various components of protostel- 
lar systems (envelope, outflow, disk). The main advantage of 
CO compared with other molecules (including water) is that its 
chemistry is simple, with most carbon locked up in CO in dense 
clouds. Also, its evaporation temperature is low, arou nd 20 K for 
pure CO ice dCollings et alj|2003l fOberg et al.ll2005l) . so that its 
freeze-out zone is much smaller than that of water. Most ground- 
based observations of CO and its isotopologues have been lim- 



ited to the lowest rotational lines originating from levels up to 
35 K. ISO has detecte d strong far-infrared C O lines up to / u =29 
from Class sources dGiannini et al.ll200ll) but the emission is 
spatially unresolved in the large 80" beam. ISO also lacked the 
spectral resolution needed to separate the shocked and quiescent 
gas or to detect intrinsically- weaker I3 CO and C I8 lines on top 
of the strong continuum. 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



T he NGC 1333 region in Perseus (d = 235 pc: lHirota et al.l 
2008) contains several deeply embedded Class s ources within 
a ~1 pc region driving pow erful outflows (e.g., ILiseau et al.l 
ll988blHatchell & FulledboOST) , The protostars IRAS 4A and 4B, 
separated by ~31", and IRAS 2 A are prominent sub millimeter 
continuum sources (luminosities of 5.8, 3.8 and 2 L Q ) with en- 
velope masses of 4.5, 2.9 and 1.0 M Q , respectively (Sandel l et al.l 
1991: Ltergensen et alj|2009l) . All three are among the brightest 
and best studied low-mass sources in terms of molecular lines , 
with several complex molecules detected (e.g jBlake et al.8 l995: 
iBottinelli etal]|2007l) . Here HIFI data of CO and its isotopo- 
logues are presented for these three sources and used to quantify 
the different physical c omponents. In an accompanying letter, 
iRristensen et al.l d2010h present complementary HIFI observa- 
tions of H2O and analyze CO/H2O abundance ratios. 
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Table 1. Overview of the observations of IRAS 2A, 4A, and 4B. Table 2. Observed line intensities. 
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f l) |j0rgensen et all <20021) ; (2) JCMT archive; (3) Yildiz et al. (in prep.); (4) this 
work 



2. Observations and results 

The N GC 1333 data were obtained with HIFI dde Graauw et al 
2010) onboard the Herschel Space Observatory dPilbratt et al 
2010), in the context of the WISH key program (van Dishoeck 



et al. in prep.). Single pointings at the three source positions were 
carried out between 2010 March 3 and 15 during the Herschel 
HIFI priority science program (PSP). Spectral lines were ob- 
served in dual-beam switch (DBS) mode in HIFI bands la, 4a, 
4b, and 5a with a chop reference position located 3' from the 
source positions. The observed positions (J2000) are: IRAS 2A: 
3 h 28 m 55 s .6, +31°14'37'.'l; IRAS 4A: ^ 29^^.5, +3ri3 / 30 , /9 ; 
and IRAS 4B: 3 h 29 m 12 s :0, +31°13'08'.'l dj0rgensen et al.ll2009l) . 

Table Q] summarizes the lines observed with HIFI together 
with complementary lower-/ lines obtained with ground-based 
telescopes. The Herschel data were taken using the wide band 
spectrometer (WBS) and high resolution spectrometer (HRS) 
backends. Owing to the higher noise ( V2 more) in HRS than 
WBS, mainly WBS data are presented here. Only the narrow 
C ls O 5—4 lines use the HRS data. Integration times (on+off) are 
10, 20, 30, 40, and 60 minutes for the 12 CO 10-9, C ls O 9-8, 
10-9, 13 CO 10-9, and C 18 5-4 lines respectively. The HIFI 
beam sizes correspond to ~20" (~4700 AU) at 1152 GHz and 
-42" (-10000 AU) at 549 GHz. Except for the 12 CO 10-9 line, 
all isotopologue lines were observed together with H2O lines. 

The calibration uncertainty for the HIFI data is of the order 
of 20% and the pointing accuracy is around 2". The measured 
line intensities were converted to the main-beam brightness tem- 
peratures 7mb = T^/t]mb by using a beam efficiency 77MB = 0.74 
for all HIFI lines. Data processing started from the standard 
HIFI pipeline in the Herschel intera ctive processing environment 
(HIPE3) ver. 3.0.1 (lOtt et al.ll20Toh . where the V LS r precision is 
of the order of a few m s _1 . Further reduction and analysis were 
done using the GILDAS-CLASjfl software. The spectra from the 
H- and V-polarizations were averaged in order to obtain a better 
S IN. In some cases a discrepancy of 30% was found between the 
two polarizations, in which case only the H band spectra were 
used for analysis since their rms is lower. 

Complementary ground-based spectral line observations 
of 12 CO 6-5 were obtained at the 12-m Atacama Pathfinder 
Experiment Te lescope (APEX), us ing the CHAMP+ 2x7 pixel 
array receiver (Gli sten et al.l l2008h . The lower-/ spectral lines 



1 HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 

2 http://www.iram.fr/IRAMFR/GILDAS/ 
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<n) In 0.5 km s" 1 bins. 



were obtained from the James Clerk Maxw ell Telescope (JCMT) 
archive and from J0rge nsen et al.l (120021) . Details will be pre- 
sented elsewhere (Yildiz et al., in prep.). 

The observed line profiles are presented in Fig.Q]and the cor- 
responding line intensities in Table|2] For the 12 CO 10-9 toward 
IRAS 2A, the emission from the blue line wing was chopped out 
due to emission at the reference position located in the blue part 
of the SVS 13 outflow. A Gaussian fitted to the red component 
of the line was used to obtain the integrated intensity. 

iRristensen et al.l (12010b identify three components in 
the H2O line profiles centered close to the source ve- 
locities: a broad underlying emission profile (Gaussian 
with FWHM ~25-30 km s~'), a medium-broad emission pro- 
file (FWHM ~5-10 km s~'), and narrow self-absorption lines 
(FWHM -2-3 km s" 1 ); see the H 2 2 02 -ln lines in Fig.Q] The 
same components are also seen in the CO line profiles, albeit less 
prominently than for H2O. The broad component dominates the 
12 CO 10-9 lines of IRAS 4A and 4B and is also apparent in the 
deep 12 CO 6-5 spectrum of IRAS 2A (Fig.0. The medium com- 
ponent is best seen in the 13 CO 10-9 profiles of IRAS 4 A and 
4B and as the red wing of the 12 CO 10-9 profile for IRAS 2A. 
A blow-up of the very high S IN spectrum of C 18 5-4 for 
IRAS 4A (insert in Fig. |T} also reveals a weak C 18 medium- 
broad profile. The narrow component is cle arly observed in 
C 18 Q emission and I2 CO low-/ self-absorption. iKristensen et al.l 
(120101) interpret the broad component as shocked gas along the 
outflow cavity walls, the medium component as smaller-scale 
shocks created by the outflow in the inner (<1000 AU) dense 
envelope, and the narrow component as the quiescent envelope, 
respectively. 
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Fig. 1. Spectra at the central positions of IRAS 2A, 4A and 4B 
Top to bottom: H2O202-I11 line from iKristensen et al.l £2010) 
illustrating the medium and broad components, and spectra of 
12 CO, 13 CO, and C ls O. The red lines correspond to the source 
velocities as obtained from the low-/ C 18 lines. The insert in 
the C ls O 5^4 line for IRAS 4A illustrates the weak medium 
component with peak Tmb = 22 mK obtained after subtracting a 
Gaussian fit to the narrow line. 



3. Analysis and discussion 

3. 1 . Broad and medium components: shocked gas 

To quantify the physical properties of the broad outflow com- 
ponent, line ratios are determined for the wings of the line pro- 
files. Figure [2] shows the CO 6-5/CO 10-9 ratio as a function 
of velocity. The APEX-CHAMP+ CO 6-5 maps of IRAS 4A 
and 4B from Yi l diz et al. (in prep.) and IRAS 2 A from 
Ivan Kempen et al.1 d2009l) are resampled to a 20" beam so that 
both lines refer to the same beam. The ratios are compared 
with model non-LT E excitation line intens ities calculated us- 
ing the RADEX code (Ivan der Tak et al.ll2007l) (Fig. \KJ\ Online 
Appendix A). The density within a 20" diam eter is taken to 
be >10 5 crrT 3 based on the modeling results of J0rg ensen et al.l 
(2002| see also Sect. 3.3 and Appen dix A). The detection of 
medium-broad CS 10-9 emission by J0r gensen et al.l d2005b) 
toward IRAS 4A and 4B indicates densities of the order of a few 
10 6 cm 3 . For the range of densities indicated in Fig. IA.1I the 
line ratios imply high temperatures: IRAS 2A, 7\m =70-130 K; 
IRAS 4A, Tkin =90-120 K; and IRAS 4B, 7\ in =140-180 K. 
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Fig. 2. Ratios of CO 6-5/CO 10-9. Top: CO line profiles. The 
CO 6-5 and 10-9 profiles have been multiplied by a factor of 
2 for IRAS 2A and 4B. Middle and bottom: ratio of line wing 
intensity in the specified velocity range indicated in the top panel 
for the red and blue wings. 



The optical depth of the 12 CO emission is constrained by the 
12 CO 10-9/ 13 CO 10-9 ratios. For IRAS 4B, the optical depth 
of the 12 CO line wings is found to drop with velocity, ranging 



from r. 



wing 



42 near the center to ~0.4 at the highest veloci- 
ties where 13 CO is detected. This justifies the assumption that 
the broad 12 CO 10-9 lines are optically thin. Total CO column 
densities in the broad component for these conditions are 4 and 
1 xlO 16 cirT 2 for IRAS 4A and 4B, respectively. For IRAS 2A, 
the broad column density is calculated from the CO 6-5 spec- 
trum as 6 XlO 15 cut 2 . Using CO/H 2 = 10~ 4 gives the H 2 col- 
umn densities listed in Table[3] 

The medium component attributed to small-scale shocks in 
the inner envelope can be probed directly by the I3 CO 10-9 data 
for IRAS 4A and 4B. For IRAS 2A, the Gaussian fit to the red 
wing of the 12 CO 10-9 is used. By assuming a similar range 
of temperatures and densities as for the broad component, beam 
averaged 12 CO column densities of 2, 6, and 2 xlO 16 cm -2 are 
found for IRAS 2A, 4A, and 4B respectively, if the lines are 
optically thin and using 12 C/ 13 C = 65. The very weak medium 
component found in the C ls O 5-4 profile for IRAS 4A agrees 
with this value if the emission arises from a compact (few ") 
source. Assuming CO/H2=10~ 4 leads to the numbers in Table[3] 
The overall uncertainty in all column densities is a factor of 2 
due to the range of physical conditions used to derive them and 
uncertainties in the adopted CO/H2 ratio and calibration. The to- 
tal amount of shocked gas is < 1 % of the total gas c olumn density 
in the beam for each source ( J0rg ensen et alj|2002l) . 

3.2. Narrow component: bulk warm envelope 

The narrow width of the C ls O emission clearly indicates an 
origin in the quiescent envelope. Naively, one would associate 
emission coming from a level with E u /ks- 237 K (9-8) with the 
warm gas in the innermost part of the envelope. To test this hy- 
pothesis, a series of envelope models was run with varying CO 
abundance profiles. The models were constructed assuming a 
power-law density structure and then calculating the temperature 
structure by fitting both the far-infrared spectral energy distribu- 
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Table 3. Summary of column densities, 7V(H2) in cm 2 in the 
broad and medium components in 20" beam. 



Source 


Broad Comp. 


Medium Comp. 


IRAS 2A 


6 X 10 1 * a 


2xl0 2Ut 


IRAS 4A 


4 x 10 20 * 


6 x 10 20 c 


IRAS 4B 


1 x 10 20 * 


2 x 10 20 c 



Obtained from (fl) CO 6-5, <b) CO 10-9, (c) B CO 10-9 spectra. 




T [K] T Q [K] 



Fig. 3. Dependence of line intensities on temperature Tq of C 18 
(left) and 13 CO (right) for an "anti-jump" model of the CO abun- 
dance in the IRAS 2A envelope. The line intensities are mea- 
sured relative to a model where the CO abundance is undepleted 
at all radii. Each curve therefore represents the fraction of the 
line intensity for the given transition, which has its origin in gas 
at temperatures below Tq. The dashed lines indicate the levels 
corresponding to 50 and 90% respectively. 

tion and the submillimeter spatial extent (J0rgensen et al. 2002). 
Figure[3]compares the fractional line intensities for the C 18 and 
13 CO transitions in a spherical envelope model for IRAS 2 A as 
a function of temperature. In these models, the abundance in the 
outer envelope was kept high, Xq=2.1 X 10~ 4 with respect to H2 
(all available gas-phase carbon in CO), decreasing by a factor of 
1000 at temperatures higher tha n a specific temper ature, Tq (a 
so-called 'anti-jump' model (see Schoier et al. 2004, for nomen- 
clature). These models thereby give an estimate of the fraction 
of the line emission for a given transition (in the respective tele- 
scope beams) which has its origin at temperatures lower than Tq. 

For C ls O, 90% of the emission in the transitions up to and 
including the 5^4 HIFI transition has its origin at temperatures 
lower than 25-30 K, meaning that these transitions are predom- 
inantly sensitive to the outer parts of the protostellar envelope. 
The 9-8 transition is more sensitive to the warm parts of the 
envelope, but still 50% of the line flux appears to come from 
the outer envelope with temperatures less than 50 K. The 13 CO 
transitions become rapidly optically thick in the outer envelopes: 
even for the 9-8 transition, 90% of the line flux can be associated 
with the envelope material with temperatures lower than 40 K. 

The C ls O 9-8 line is clearly a much more sensitive probe 
of a CO ice evaporation zone than any other observed CO line. 
J0rg ensen et al.l d2005c) showed that the low-/ C ls O lines re- 
quire a drop in the abundance at densities higher than 7 x 10 4 
cm -3 due to freeze-out. However, they did not have strong proof 
for CO evaporation in the inner part from that dataset. Using 
the temperature and density structure for IRAS 2A as described 
above, we computed the C I8 line intensities in the respec- 
tive tele scope beams following the method by J0rg ensen et al.l 
(2005c). In this 'anti-jump' model, the outer C 18 abundance 



is kept fixed at Xq = 5.0 x 10~ 7 , whereas the inner abundance 
Xq and the freeze-out density rife are free parameters. A x 1 fit 
to only the C ls O 1-0, 2-1 and 3-2 lines gives best-fit values of 
X n — 3 X 10~ 8 and = 7 x 10 4 cm 3 , consistent with those 
of J0rgen sen et al.l d2005d) . However, this model underproduces 
the higher-/ lines by a factor of 3^1 (Fig. IB.2l in Appendix B). 

To solve this underproduction, the inner abundance has to 
be increased in a so-called 'drop-abundance' profile. The fit pa- 
rameters are now the inner abundance X m and the evaporation 
temperature T ev , keeping Xq and «de fixed at the above values. 
Figure lB~5l in Appendix B shows the^ 2 plots to the C 18 6-5 and 
9-8 lines. The evaporation temperature is not well constrained, 
but low temperatures of T ev « 25 K are favored because they pro- 
duce more C ls O 5^4 emission. The best-fit X in - 1.5 x 10~ 7 in- 
dicates a jump of a factor of 5 compared with Xu- Alternatively, 
T ev can be kept fixed at 25 K and both X m and Xu can be var- 
ied by fitting all five lines simultaneously. In this case, the same 
best-fit value for Xi n is found but only an upper limit on Xu of 
~ 4 x 10~ 8 . Thus, for this physical model, X{ n > Xu, implying 
that a jump in the abundance is needed for IRAS 2A. 

4. Conclusions 

Spectrally resolved Herschel-HIFI observations of high-/ CO 
lines up to 12 CO 10-9 and C ls O 9-8 have been performed to- 
ward three low-mass young stellar objects for the first time. 
These data provide strong constraints on the density and tem- 
perature in the various physical components, such as the quies- 
cent envelope, extended outflowing gas, and small-scale shocks 
in the inner envelope. The derived column densities and tem- 
peratures are important for comparison with water and other 
molecules such as O2, for which HIFI observations are planned. 
Furthermore, it is shown conclusively that in order to reproduce 
higher-/ C 18 lines within the context of the adopted physical 
model, a jump in the CO abundance due to evaporation is re- 
quired in the inner envelope, something that was inferred, but not 
measured, from ground-based observations. Combination with 
even higher-/ CO lines to be obtained with Herschel-PACS in 
the frame of the WISH key program will allow further quantifi- 
cation of the different physical processes invoked to explain the 
origin of the high-/ emission. 
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Table B.l. Summary of CO abundance profiles for IRAS 2A. 
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Fig. A.l. Model line ratios of CO 6-5/10-9 for a slab model with 
a range of temperatures and densities. The adopted CO column 
density is 10 17 cuT 2 with a line width of 10 km s _1 , comparable 
to the inferred values. For these parameters the lines involved 
are optically thin. The colored lines give the range of densities 
w ithin the 20" beam for the three sources based on the models 
of lJ0rgensen et al.1 d2002l) . 



Appendix A: Radex model 



Figure |AT| shows the CO 6-5/10-9 line ratios for a slab model 
with a range of temperatures and densities. 



Appendix B: Abundance profiles for IRAS 2A 

Among the three sources, IRAS 2A has been selected for de- 
tailed CO abundance profile modeling because more data are 
available on this source, and because its physical and chem- 
ical structure has been well characterized through the high 
angular resolution submillimeter single dish and interferomet- 
ric observations of iJOrgensen etail 00021 l2005al) . The physi- 
cal param^ters_are taken from the continuum modeling results 
of lj0rgensen et al.1 d2002). In that paper, the ID dust radiative 
transfer code DUSTY rilvezic&Elitziui ri997) was used assum- 
ing a power law to describe the density gradient. The dust tem- 
perature as function of radius was calculated self-consistently 
through radiative transfer given a central source luminosity. 
Best-fit model parameters were obtained by comparison with 
the spectral energy distribution and the submillimeter contin- 
uum spatial extent. The resulting envelope structure parameters 
are used as input to the Ratran rad iative transfer modeling code 
(Hogerheiide & van der Tak 2000) to model the CO line inten- 
sities for a given CO abundance structure through the envelope. 
The model extends to 11000 AU from the protostar, where the 
density has dro pped to 2 X IP 4 cm" 3 . The CO-H2 collisional rate 
coefficients of Yang et al] d2010h have been adopted. 

The C 18 lines are used to determine the CO abundance 
structure because the lines of this isotopologue are largely op- 
tically thin and because they have well-defined Gaussian line 
shapes originating from the quiescent envelope without strong 
contaminations from outflows. Three types of abundance pro- 
files are examined, namely 'constant', 'anti-jump' and 'drop' 
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abundance profiles. Illustrative models are shown in Fig. IB. H and 
the results from these models are summarized in Table IB. ll 



10 



10 



io-° 

v 10" 6 
o 

a 7 

3 ID" 7 
a 

5 to" 8 

< 10" 6 



10 
10" 



Constant Abundance Profile 



Anti-Jump Profile 



Drop Abundance Profile 



10' 



10 r , 

Radius [AU] 



10 4 



Fig. B.l. Examples of constant, anti-jump, and drop abundance profiles 
for IRAS 2A for T ev =25 K and « dc =7x 10 4 cnT 3 



B.1. Constant abundance model 

The simplest approach is to adopt a constant abundance across 
the entire envelope. However, with this approach, and within the 
framework of the adopted source model, it is not possible to 
simultane ously reproduce all line intensities. This was already 
shown by lj0rgensen et al.l d2005ch . For lower abundances it is 
possible to reproduce the lower-/ lines, while higher abundances 
are required for higher-/ lines. In Fig. IB.2I the C I8 spectra 
of a constant-abundance profile are shown for an abundance of 
X =1.4xl0~ 7 , together with the observed spectra of IRAS 2A. 
Based on these results, the constant-abundance profile is ruled 
out for all three sources. 



B.2. Anti-jump abundance models 

The anti-jump model is commonly adopted in mod- 
els o f pre-stellar cores w it hout a cent r al he ating source 
(e.g.. iBergin & Snelll " 



2002; iTafalla et all I2004 . Following 



J0rgensen et al.l (12005c). an anti-jump abundance profile was 
employed by varying the desorption density, n<je> an d inner abun- 
dance X{ n =Xo in order to find a fit to our observed lines. Here, 
the outer abundance Xq was kept high at 5.0xl0~ 7 correspond- 
ing to a 12 CO abundance of 2.4xl0 -4 for 16 O/ 1 8 O=550 as was 
found a ppropriate for the case of IRAS 2A by J0rg ensen et al.l 
(l2005d) . This valu e is consistent w ith the CO/H2 abundance 
ratio determined bv lLacv etal1 dl994) for dense gas without CO 
freeze-out. 

The best fit to the three lowe st C 18 Q lines (1-0, 2-1 a nd 3-2) 
is consistent with that found by J0rge nsen et al.l (12005 d) . corre- 
sponding to «d e =7xl0 4 cnT 3 and Xd=3x10 -8 (CO abundance 
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Fig. B.2. Best fit constant (green), anti-jump (blue) and drop abun- 
dance (red) Ratran models overplotted on the observed spectra. All 
spectra refer to single pointing observations. The calibration uncertainty 
for each spectrum is around 20-30% and is taken into account in the x 2 
fit. See Table lB. il for parameters. 



of 1.7 x 10~ 5 ). In the;f 2 fits, the calibration uncertainty of each 
line (ranging from 20 to 30%) is taken into account. These mod- 
eled spectra are overplotted on the observed spectra in Fig. IB. 21 
as the blue lines, and show that the anti-jump profile fits well the 
lower-/ lines but very much underproduces the higher-/ lines. 

The value of Xq was verified a posteriori by keeping «de at 
two different values of 3.4xl0 4 and 7xl0 4 cm 4 . This is illus- 
trated in Fig. IB. 3 1 where the^ 2 contours show that for both values 
of Wde, the best-fit v alue of Xq is ~5xl0~ 7 , the value also found 
in |j0rgensen et alJ d2005cl) . The x 2 contours have been calcu- 
lated from the lower-/ lines only, as these are paramount in con- 
straining the value of Xq. Different^ 2 plots were made, where it 
was clear that higher-/ lines only constrain X^, as expected. The 
effect of «de is illustrated in Fig. IB .41 for the two values given 
above. 



B.3. Drop-abundance profile 

In order to fit the higher-/ lines, it is necessary to employ a drop- 
abundance structure in which the inner abundan ce X m increases 
above the ice evaporation temperature T ev dj0rgensen et alJ 
2005c). The abundances Xd and Xq for T < T ev are kept the 
same as in the anti-jump model, but X m is not necessarily the 
same as Xq. In order to find the best-fit parameters for the higher- 
/ lines, the inner abundance X m and the evaporation temperature 
T ev were varied. The^ 2 plots (Fig. IB.5l left panel) show best-fit 
values for an inner abundance of X[ n = 1.5xl0~ 7 and an evapo- 
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Fig. B.3. The x 1 plots for the anti-jump profiles where Xo and Xn val- 
ues are varied. Right: for n de =7x! 4 and left: for n de =3.4 xl0 4 cm 3 . 
The asterisk indicates the value for Jetrgen sen et al.l l[2P05c) used here. 
Contours are plotted at the 2cr, 3<r, and 4<x confidence levels (left) and 
3<x and 4<r confidence levels (right). 
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Fig. B.4. The IRAS 2A s pectra for the X and X parameters corre- 
sponding to the values in J0rge nsen et alJ i2005d) for different /i de val- 
ues of 3.4xl0 4 and7xl0 4 cnr 3 . 
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Fig. B.5. Reduced^ 2 plots and best-fit parameters (indicated with *) 
for the anti-jump model fit to the lines of C ls O 1-0, 2-1, 3-2, 6-5 and 
9-8 (right) and for the drop abundance model fit to the higher-./ lines 
of C l8 6-5 and 9-8 (left). Contours are plotted at the Icr, 2a, 3cr, and 
4<r confidence levels. 



ration temperature of 25 K (consistent with the laboratory val- 
ues), although the latter value is not strongly constrained. These 
parameters fit well the higher-/ C I8 6-5 and 9-8 lines (Fig. 
IB.2I >. The C ls O 5-4 line is underproduced in all models, likely 
because the larger HIFI beam picks up extended emission from 
additional dense material to the northe ast of the source seen in 
BIMA C ls O 1-0 (Volgen au etalll2006h map. 
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Because the results do not depend strongly on T ev , an alter- 
native approach is to keep the evaporation temperature fixed at 
25 K and vary both X m and Xo by fitting both low- and high-/ 
lines simultaneously. In this case, only an upper limit on Xo of 
~ 4 x 1CT 8 is found (~Fig. IB.5l right panel), whereas the inferred 
value of X m is the same. This figure conclusively illustrates that 
^in > ^d, i-e-, that a jump in the abundance due to evaporation 
is needed. 

The above conclusion is robust within the context of the 
adopted physical model. Alternatively, one coul d inves tigate 
different physical models such as those used by IChiang et alj 
(2008), which have a density enhancement in the inner envelope 
due to a magnetic shock wall. This density increase could partly 
mitigate the need for the abundance enhancement although it is 
unlikely that the density jump is large enough to fully compen- 
sate. Such models are outside the scope of this paper. An ob- 
servational test of our model would be to image the N2H + 1-0 
line at high angular resolution: its emission should drop in the 
inner ~900 AU (~4") where N2H + would be destroyed by the 
enhanced gas-phase CO. 



